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We show how generalised unitarity cuts in D = 4 — 2e dimensions can be used to calcu- 
late efficiently complete one-loop scattering amplitudes in non- super symmetric Yang-Mills 
theory. This approach naturally generates the rational terms in the amplitudes, as well 
as the cut-constructible parts. We test the validity of our method by re-deriving the one- 
loop ++++, — h++, h+, — I h and gluon scattering amplitudes using 

generalised quadruple cuts and triple cuts in D dimensions. 
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1 Introduction 



Over the past year, major progress in the calculation of scattering amplitudes in per- 
turbative Yang-Mills theory has been made. This was triggered by Witten's discovery 
that tree-level amplitudes in Yang-Mills can equivalently be derived via a string theory 
calculation, where the string theory in question is the topological B model with target 
space a supersymmetric version of Penrose's twistor space [1]. Witten also observed that 
tree-level scattering amplitudes, when Fourier transformed to twistor space, have an inter- 
esting geometrical structure, namely they have support on algebraic curves; for the simple 
case of the maximally helicity violating (MHV) amplitude, described by the Parke- Taylor 
formula, the curve is just a line (for real twistor space). This remarkable observation 
gives an explanation for the unexpected and previously rather mysterious simplicity of 
tree-level scattering amplitudes in Yang-Mills such as the Parke- Taylor formula, which is 
not at all apparent in a calculation performed using standard Feynman rules. 

On a different line of development, the simplicity of tree-level scattering amplitudes 
was linked to the existence of novel recursion relations discovered by Britto, Cachazo 
and Feng (BCF) [2], and subsequently proved by the same authors and Witten (BCFW) 
[3]. The elegant proof of [3] is based on very general properties of amplitudes, such 
as analyticity [4-6] and factorisation on multiparticle poles, and hence gave rise to the 
hope that recursion relations may arise in very different contexts. Indeed, novel recursion 
relations were also found in general relativity [7, 8] , scalar theory [7] , for the finite rational 
amplitudes at one-loop in Yang-Mills and massless QCD [9,10], and for tree amplitudes 
involving massive scalars and gluons in Yang-Mills [11]. 

The simplicity of tree-level amplitudes in Yang-Mills was exploited by Bern, Dixon, 
Dunbar and Kosower (BDDK) in order to build one-loop scattering amplitudes [12,13]. 
By applying unitarity at the level of amplitudes, rather than Feynman diagrams, these 
authors were able to construct many one-loop amplitudes in supersymmetric theories, 
such as the infinite sequence of MHV amplitudes in M = A and in M = 1 super Yang-Mills 
(SYM). The unitarity method of BDDK by-passes the use of Feynman diagrams and its 
related complications, and generates results of an unexpectedly simple form; for instance, 
the one-loop MHV amplitude in M = 4 SYM is simply given by the tree-level expression 
multiplied by a sum of "two-mass easy" box functions, all with coefficient one. As a side 
remark, we would like to mention that higher-loop amplitudes in A/" =4 SYM also display 
intriguing regularities [14-16]. 

The geometrical structure in twistor space of the amplitudes was also the root of a 
further important development. In [17], Cachazo, Svrcek and Witten (CSW) proposed 
a novel perturbative expansion for on-shell amplitudes in Yang-Mills, where the MHV 
amplitudes are lifted to vertices, joined by simple scalar propagators in order to form 
amplitudes with an increasing number of negative helicities. Applications at tree level 
confirmed the validity of the method and led to the derivation of various new amplitudes 
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in gauge theory [17-24]. 



In [17], a heuristic derivation of the CSW method was given from the twistor string 
theory. Rather unfortunately, the latter only appears to describe the scattering amplitudes 
of Yang-Mills at tree level [25], as at one loop states of conformal supergravity enter the 
game, and cannot be decoupled in any known limit. The duality between gauge theory and 
twistor string theory is thus spoiled by quantum corrections. Surprisingly, it was found 
by three of the present authors that the MHV method at one-loop level does succeed in 
correctly reproducing the scattering amplitudes of the gauge theory [26]. Furthermore, 
the twistor space picture of one-loop amplitudes is now in complete agreement with that 
emerging from the MHV methods, which suggests that the amplitudes at one loop have 
localisation properties on unions of lines in twistor space; an initial puzzle [27] was indeed 
clarified and explained in terms of a certain "holomorphic anomaly", introduced in [28], 
and further analysed in [29-33]. A proof of the MHV method at tree level was finally 
given in [3]; at loop level, however, it remains a (well-supported) conjecture. 

The initial successful application of the MHV method to J\f=4 SYM [26] was followed 
by calculations of MHV amplitudes in M = 1 SYM [34, 35], and in pure Yang-Mills 
[36], where the four-dimensional cut-constructible part of the infinite sequence of MHV 
amplitudes was derived. However, amplitudes in non-supersymmetric Yang-Mills theory 
also have rational terms which escape analyses based on MHV diagrams at one loop [26,36] 
or four- dimensional unitarity [12,13]. 

Amplitudes in supersymmetric theories are of course special. They do contain rational 
terms, but these are uniquely linked to terms which have cuts in four dimensions. In 
other words, these amplitudes can be reconstructed uniquely from their cuts in four- 
dimensions [12, 13] - a remarkable result. These cuts are of course four- dimensional tree- 
level amplitudes, whose simplicity is instrumental in allowing the derivation of analytic, 
closed-form expressions for the one-loop amplitudes. In non-supersymmetric theories, 
amplitudes can still be reconstructed from their cuts, but on the condition of working in 
4 — 2e dimensions, with e =^ [37-39]. This is a powerful statement, but it also implies 
the rather unpleasant fact that one should in principle work with tree-level amplitudes 
involving gluons continued to 4 — 2e dimensions, which are not simple. 

An important simplification is offered by the well-known supersymmetric decomposi- 
tion of one-loop amplitudes of gluons in pure Yang-Mills. Given a one-loop amplitude A g 
with gluons running the loop, one can re-cast it as 

Ag = (A s + 4A f + 3A S ) - 4(A f + A s ) + A s . (1.1) 

Here At (A B ) is the amplitude with the same external particles as A g but with a Weyl 
fermion (complex scalar) in the adjoint of the gauge group running in the loop. This 
decomposition is useful because the first two terms on the right hand side of (|l.lj) are 
contributions coming from an J\f = 4 multiplet and (minus four times) a chiral Af= 1 mul- 
tiplet, respectively; therefore, these terms are four- dimensional cut-constructible, which 
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simplifies their calculation enormously. The last term in (jl.lj) . A s , is the contribution 
coming from a scalar running in the loop. The key point here is that the calculation of 
this term is much easier than that of the original amplitude A g . It is this last contribution 
which is the focus of this paper. 

The root of the simplification lies in the fact that a massless scalar in 4 — 2e dimensions 
can equivalently be described as a massive scalar in four dimensions [38,39]. Indeed, if L 
is the (4 — 2e)-dimensional momentum of the massless scalar (L 2 = 0), decomposed into a 
four-dimensional component Z(4) and a — 2e-dimensional component l(-2e), L := l(4)+l(-2e), 
one has L 2 := l 2 ^ + lf_ 2e ) — I a) ~ ^ \ where l 2 _ 2e ) := an d the four- dimensional and 
— 2 e- dimensional subspaces are taken to be orthogonal. The tree-level amplitudes entering 
the (4— 2e)-dimensional cuts of a one- loop amplitude with a scalar in the loop are therefore 
those involving a pair of massive scalars and gluons. Crucially, these amplitudes have a 
rather simple form. Some of these amplitudes appear in [38,39]; furthermore, a recent 
paper [11] describes how to efficiently derive such amplitudes using a recursion relation 
similar to that of BCFW. 

Using two-particle cuts in 4 — 2e dimensions, together with the supersymmetric decom- 
position mentioned above, various amplitudes in pure Yang-Mills were derived in recent 
years, starting with the pioneering works [38,39]. In this paper we show that this analysis 
can be performed with the help of an additional tool: generalised (4 — 2e)-dimensional 
unitarity. 

Generalised four-dimensional unitarity [5,6,40-42] was very efficiently applied in [43] 
to the calculation of one- loop amplitudes in Af = 4 SYM. Amplitudes in this theory can be 
written as a sum of box functions, multiplied by rational coefficients. To each box function 
is uniquely associated a (generalised) quadruple cut, so that, schematically, each coefficient 
of a box function is expressed as a particular quadruple cut of the one-loop amplitude, 
which is nothing but a product of four tree-level amplitudes. Generalised cuts require 
the amplitudes to be continued to complexified Minkowski space, which in turn has the 
consequence that three-point amplitudes no longer vanish, and enter the cut-amplitude in 
an important way [43]. 1 The calculation of one-loop amplitudes in M =4 SYM was in this 
way turned into an algebraic problem [43] . Using generalised unitarity in four dimensions, 
the infinite sequence of next-to-MHV amplitudes in M = 4 SYM was determined [44] ; 
generalised unitarity was also applied to M = 1 SYM, in particular to the calculation of 
the next-to-MHV amplitude with adjacent negative-helicity gluons [45]. These amplitudes 
can be expressed solely in terms of triangles, and were efficiently computed in [45] using 
triple cuts. 2 

The main point of this paper is the observation that generalised unitarity is actually 
a useful concept also in 4 — 2e dimensions; in turn this means that generalised (4 — 2e)- 

^^This circumstance extends to the (4 — 2e)-dimensional three-point scattering amplitudes which will 
be considered in this paper. 

2 A new calculation based on localisation in spinor space was also introduced in [46]. 
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dimensional unitarity is relevant for the calculation of non-supersymmetric amplitudes 
at one loop. In particular in this paper we will be able to compute amplitudes in non- 
supersymmetric Yang-Mills by using quadruple and triple cuts in 4 — 2e dimensions. 
This is advantageous for at least three reasons. First of all, working with multiple cuts 
simplifies considerably the algebra, because several on-shell conditions can be used at 
the same time; furthermore, for the case of quadruple cuts the integration is actually 
completely frozen [43] so that the coefficient of the relevant box functions entering the 
amplitude can be calculated without performing any integration at all. Lastly, the tree- 
level sub- amplitudes which are sewn together in order to form the multiple cut of the 
amplitude are simpler than those entering the two-particle cuts of the same amplitude. 
It seems clear that immediate further progress with this approach will not require major 
new conceptual advances, and that it will be directly applicable to more complicated and 
currently unknown amplitudes. 

We describe this method in some detail in Section 2, and then move on to present 
various examples of its application. Specifically, using generalised unitarity in 4 — 2e 
dimensions we will re-calculate the all-orders in e expressions of all one-loop, four gluon 
scattering amplitudes in non-supersymmetric Yang-Mills, that is ++++, — and 

the two MHV amplitudes h+ and — I h; and finally, the five-gluon all-plus helicity 

amplitude +++++. These amplitudes have already been computed to all orders in e 
in [38], and we find in all cases complete agreement with the results of that paper. The 
examples we consider are complementary, as they show that this method can be applied to 
finite amplitudes without infrared divergences, as well as to infrared divergent amplitudes 
containing both rational and cut-constructible terms. These calculations are described in 
Section 3 and Section 4. In an Appendix we have collected some useful definitions and 
formulae. 



2 Generalised Unitarity in D = 4 — 2e Dimensions 



Conventional unitarity and generalised unitarity in four dimensions have been shown to be 
extremely powerful tools for calculating one-loop and higher-loop scattering amplitudes in 
supersymmetric gauge theories and gravity. At one-loop, conventional unitarity amounts 
to reconstructing the full amplitude from the knowledge of the discontinuity or imaginary 
part of the amplitude. In this process the amplitude is cut into two tree-level, on-shell 
amplitudes defined in four dimensions, and the two propagators connecting the two sub- 
amplitudes are replaced by on-shell delta-functions which reduce the loop integration to 
a phase space integration. In principle this cutting technique is only sensitive to terms 
in the amplitude that have discontinuities, like logarithms and polylogarithms, and in 
general any cut-free, rational terms are lost. However, in supersymmetric theories all 
rational terms turn out to be uniquely linked to terms with discontinuities, and therefore 
the full amplitudes can be reconstructed in this fashion [12, 13]. 
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Furthermore, in supersymmetric theories the one-loop amplitudes are known to be 
linear combinations of scalar box functions, linear triangle functions and linear bubble 
functions, with the coefficients being rational functions in spinor products. So the task is 
really to find an efficient way to fix those coefficients with as few manipulations and/or 
integrations as possible. 

The method based on conventional unitarity introduced by BDDK in [12, 13] does 
not evaluate the phase space integrals explicitly (from which the full amplitude would be 
obtained by performing a dispersion integral), rather it reconstructs the loop integrand 
from which one is able to read off the coefficients of the various integral functions. In 
practice this means that for a given momentum channel the integrand (which is a product 
of two tree amplitudes) is simplified as much as possible using the condition that the two 
internal lines are on-shell, and only in the last step the two delta-functions are replaced by 
the appropriate propagators which turn the integral from a phase space integral back to a 
fully-fledged loop integral. The resulting integral function will have the correct disconti- 
nuities in the particular channel, but, in general, it will also have additional discontinuities 
in other channels. Nevertheless, working channel by channel one can extract linear equa- 
tions for the coefficients which allow us in the end to determine the complete amplitude. 
However, because of the problem of the additional, unwanted discontinuities, this does 
not provide a diagrammatic method, i.e. one cannot just sum the various integrals for 
each channel since different discontinuities might be counted with different weights. 

It is natural to contemplate if there exist other complementary, or more efficient meth- 
ods to extract the above mentioned rational coefficients of the various integral functions, 
and if in particular we can replace more than two propagators by delta functions, so that 
the loop integration is further restricted - or even completely localised. The procedure of 
replacing several internal propagators by 5^ + )-functions is well known from the study of 
singularities and discontinuities of Feynman integrals, and goes under the name of gen- 
eralised unitarity [5,6]. What turns generalised unitarity into a powerful tool is the fact 
that generalised cuts of amplitudes can be evaluated with less effort than conventional 
two-particle cuts. 

The most dramatic simplification arises from using quadruple cuts in one-loop ampli- 
tudes in M = 4 SYM. In this case it is known that the one-loop amplitudes are simply 
given by a sum of scalar box functions without triangles or bubbles [12]. Each quadruple 
cut singles out a unique box function, and because of the presence of the four ^W-functions 
the loop integration is completely frozen; hence, the coefficient of this particular box is 
simply given by the product of four tree-level scattering amplitudes [43] . An important 
subtlety arises here because quadruple cuts do not have solutions in real Minkowski space; 
therefore at intermediate steps one has to work with complexified momenta. 

At this point we can push the analogy with the "reconstruction of the Feynman inte- 
grand" a step further. Using the on-shell conditions we can pull out the prefactor which 
is just the product of four tree-level amplitudes in front of the integral, and the integrand 
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of the remaining loop integral becomes just a product of four <5^ + ^-functions. If we now 
promote the integral to a Feynman integral by replacing all 5( + )-functions by the cor- 
responding propagators 3 we arrive at the integral representation of the appropriate box 
function. Note that no overcounting issue arises, because each quadruple cut selects a 
unique box function, and the final result is obtained by summing over all quadruple cuts. 
In some sense, one can really think of this as a true diagrammatic prescription. 

As we reduce the amount of supersymmetry to M— 1, life becomes a bit more compli- 
cated, since the one-loop amplitudes are linear combinations of scalar box, triangle and 
bubble integral functions. No ambiguities related to rational terms occur however, thanks 
to supersymmetry. It is therefore natural to attack the problem in two steps: First, use 
quadruple cuts to fix all the box coefficients as described in the previous paragraph. Sec- 
ond, use triple cuts to fix triangle and bubble coefficients. Note that the triple cuts also 
have contributions from the box functions which have been determined in the first step. 
The three -functions are not sufficient to freeze the loop integration completely, and 
it is advantageous to use again the "reconstruction of the Feynman integrand" method, 
i.e. use the on-shell conditions to simplify the integrand as much as possible, and lift the 
integral to a full loop integral by reinstating three propagators. The resulting integrand 
can be written as a sum of (integrands of) scalar boxes, triangles and bubbles, after 
standard reduction techniques, like Passarino-Veltman, have been employed. 

At this point it is useful to distinguish three types of triple cuts according to the 
number of external lines attached to each of the three tree-level amplitudes. If p of the 
three amplitudes have more than one external line attached, we call the cut a p-mass 
triple cut. Let us start with the 3- mass triple cut. The box terms can be dropped as 
they have been determined using quadruple cuts, the coefficients of three-mass triangles 
can be read off directly, and the remaining terms, which are bubbles or triangles with a 
different triple cut, are dropped as well. Special care is needed for 1-mass and 2- mass 
triple cuts. First let us note that any bubble can be written as a linear combination 
of scalar and linear 1-mass triangles or scalar and linear 2-mass triangles depending on 
whether the bubble depends on a two-particle invariant, if = (pi + Pi+i) 2 , or on a r- 
particle invariant, tf^ = (pi + . . . + Pi+ r -i) 2 , with r > 2. Therefore, what we want to 
argue is that two-particle cuts are not needed and that 1-mass, 2-mass and bubbles can 
be determined from the 1-mass and 2-mass triple cuts. Now every 1-mass triple cut is in 
one-to-one correspondence with a unique two-particle channel tf^ = {pi+Pi+i) 2 and allows 
us to extract the coefficients of 1-mass triangles and bubbles by only keeping terms in the 
integral depending on that particular tj 2 ' and dropping all boxes and triangles/bubbles 
not depending on that particular variable. The 2-mass triple cut is associated with two 
momentum invariants, say P 2 and Q 2 , and we only keep 2-mass triangles and bubbles 
that depend on those two invariants. 

In non-supersymmetric theories we have to face the problem that the amplitudes con- 
3 We thank David Kosower for discussions on this point. 



6 



tain additional rational terms that are not linked to terms with discontinuities. This 
statement is true if we only keep terms in the amplitude up to O(e ). If we work however 
in D = 4 — 2e dimensions and keep higher orders in e, even rational terms R develop dis- 
continuities of the form R(— s)~ e = R — elog(—s)R+0(e 2 ) and become cut-constructible 4 . 
In practice, this means that, in our procedure, whenever we cut internal lines by replacing 
propagators by <5^ + ^-functions we have to keep the cut lines in D dimensions, and in order 
to proceed we need to know tree-amplitudes with two legs continued to D dimensions. 
Because of the supersymmetric decomposition of one-loop amplitudes in pure Yang-Mills, 
which was reviewed in the Introduction, we only need to consider the case of a scalar 
running in the loop. Furthermore, the massless scalar in D dimensions can be thought 
of as a massive scalar in four dimensions L 2 = l 2 ^ + l 2 _ 2t ) — ^ 2 4 )— whose mass has to 
be integrated over [38,39]. Interestingly, a term in the loop integral with the insertion of 
"mass" term (/i 2 )" 1 can be mapped to a higher-dimensional loop integral in 4 + 2m — It 
dimensions with a massless scalar [38,39]. Some of the required tree amplitudes with 
two massive scalars and all positive helicity gluons have been calculated in [38, 39] using 
Feynman diagrams and recursive techniques, and more recently all amplitudes with up to 
four arbitrary helicity gluons and two massive scalars have been presented in [11]. 

The comments in the last paragraph make it clear that generalised unitarity techniques 
can readily be generalised to D dimensions and be used to obtain complete amplitudes 
in pure Yang-Mills and, more generally, in massless, non- supersymmetric gauge theo- 
ries. The integrands produced by the method described for four dimensional unitarity 
will now contain terms multiplied by {fi 2 ) m and, therefore, the set of integral functions 
appearing in the amplitudes includes, in addition to the four-dimensional functions, also 
higher-dimensional box, triangle and bubble functions (some explicit examples of higher- 
dimensional integral functions can be found in Appendix A). For example the one-loop 
++++ gluon amplitude, which vanishes in SYM, is given by a rational function times 
a box integral with /i 4 inserted, h[fj> 4 } = (— e)(l — e)lf~ 2€ = —1/6 + 0(e). Hence this 
amplitude is a purely rational function in spinor variables. 

In the following sections we will describe in detail how this procedure is applied in 
practice by recalculating all four-gluon scattering amplitudes and the positive helicity five- 
gluon scattering amplitude in pure Yang-Mills at one-loop level. These examples include 
the cases of infrared finite amplitudes that are purely rational (and their supersymmetric 
counterparts vanish), and infrared divergent amplitudes that contain both rational and 
cut-constructible terms. 

4 The idea of using unitarity in D = 4 — 2e dimensions goes back to [37], and was used in [38, 39]. 
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3 Four-point amplitudes in pure Yang-Mills 



In this section we recalculate all the known four-gluon scattering amplitudes, that is 
++++, — h++, h+, and finally — I h, from quadruple and triple cuts. 



3.1 The one- loop ++++ amplitude 



The one-loop ++++ amplitude with a complex scalar running in the loop is the simplest 
of the all-plus gluon amplitudes, and was first derived in [47] using the string-inspired 
formalism. 

The expression in 4 — 2e dimensions, valid to all-orders in e, is computed in [38] and 
is given by 

A scalar / 1 + 9+ q+ A+\ 2 * I 12 ]! 34 ] 

4 (1 ' 2 ' 3 ' 4 } " (470^(12X34) Ka ' (3 ' 1} 

where 5 

K, := J 4 [/i 4 ] = - c (l - c)jf*-* = -\ + 0(e). (3.2) 

o 

In this paper we closely follow the conventions of [38], with 
^ =4 " 2e [/(P,/i 2 )] := ^(-) n+1 (47r) 2 - f dH fil '^ 



(27T)4 (27T)-2e ( p _ fx 2 ) ... [{l _ £-1 Ki) 2 _ ^2] ' 

(3.3) 

where K,i are external momenta (which, in colour-ordered amplitudes, are sums of ad- 
jacent null momenta of the external gluons) and f(l,fi 2 ) is a generic function of the 
four-dimensional loop momentum I and of fi 2 . 

The amplitude with four positive helicity gluons is part of the infinite sequence of 
all-plus helicity gluons, for which a closed expression was conjectured in [48,49]. The 
result for all n is given by 

a ,\ _±_ sr (ia\ 

nl+ ''"' +j ' 48vr 2 L> (12)<23)-..(nl) ' 1 1 

or, alternatively, 

a i X "> ^1112^1314 ^1113^214 ^1114^12*3 4z€(iii2^3^4) /q p.\ 

" 1 ~96V 2^ (12) (23) ■■■(nl) ' 1 j 

I<i 1 <i2<h<i4,<n N ' 

where e(abcd) := e flups a' J 'b l/ c p d a . As e — > 0, (|3.1|) becomes 

4 48tt 2 (12) (23) (34) (41) ' 1 J 



5 Notice also that [12][34]/((12)(34)) = -si 2 s 23 /((12)(23)(34)(41)). 
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We see that this amplitude (j3.1j) consists of purely rational terms, which are cut-free in 
four dimensions. We now show how to derive (j3.1|) from quadruple cuts in D = 4 — 2e 
dimensions. 



2 




u 




-a- 





Figure 1: One of the two quadruple- cut diagrams for the amplitude l + 2 + 3 + 4 + . This 
diagrams is obtained by sewing tree amplitudes (represented by the blue bubbles) with an 
external positive-helicity gluon and two internal scalars of opposite "helicities" . There 
are two such diagrams, which are obtained one from the other by flipping all the internal 
helicities. These diagrams are equal so that the full result is obtained by doubling the 
contribution from the diagram in this Figure. The same remark applies to all the other 
diagrams considered in this paper. 

Consider the quadruple-cut diagram in Figure 1, which is obtained by sewing four 
three-point scattering amplitudes 6 with one massless gluon and two massive scalars of 
mass /i 2 . From [11] we take the three-point amplitudes for one positive-helicity gluon and 
two scalars: 



where l\ + li + k = 0. Here \q) is an arbitrary reference spinor not proportional to \k). It 
is easy to see [11] that ()3.7|) is actually independent of the choice of \q). 

The Z)-dimensional quadruple cut of the amplitude +++-1- is obtained by combining 
four three-point tree-level amplitudes, 



6 In the following for the purpose of calculating the (generalised) cuts we drop factors of i appearing 
in the usual definition of tree-amplitudes and propagators. For quadruple and two-particle cuts this does 
not affect the final result, while for triple cuts this introduces an extra (—1) factor which we reinstate at 
the end of every calculation. 




(3.7) 



(giNl] (ggNj] <g 3 |Z 3 |3] (g 4 |k|4] 
(gil) (? 2 2> (<fe3) (g 4 4) 



(3.8) 
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and similarly 



Figure 2: One of the possible three-particle cut diagrams for the amplitude l + 2 + 3 + 4 + . 
The others are obtained from this one by cyclic relabeling of the external particles. 

The reference momenta i — 1, . . . , 4 in each of the four ratios in this expression may 
be chosen arbitrarily. Then, using momentum conservation, 

l 2 — l 1 -k 2 , h = h- k A , (3.9) 

the fact that the external momenta are null, and that the internal momenta square to /i 2 , 
it is easy to see that 

ijtm = V H , (3 . 10) 

<?il> M * (12) ' 1 J 

(ft3> (g 4 4) " " (34) ' (3 ' n) 
so that the above expression ()3.8|) becomes simply 

[12][34 J (3 12) 

M (12)(34) • 

Finally, we lift the quadruple-cut box to a box function by reinstating the appropriate 
Feynman propagators. These propagators then combine with the additional factor of fi A 
in (J3.12j) to yield the factor iK±/ ^47r) 2 ~ e which is proportional to the scalar box integral 
defined in ()3.2j) . Including an additional factor of 2 due to the fact that there is a complex 
scalar propagating in the loop, the amplitude (|3.1|) is reproduced correctly. 

Next we inspect three-particle cuts. One of the three tree-level amplitudes we sew 
in the triple-cut amplitude is an amplitude with two posit ive-helicity gluons and two 
scalars [39] 

^l(/+,l + ,2 + ,/ 2 -) = V 2 {imh + ki)2 _^ 2] ■ (3-13) 
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Consider, for example, the three-particle cut defined by 1 + , 2 + , (3 + , 4 + ), see Figure 2. 
Using ()3.7|) and (|3.13|) . the product of the three tree-level amplitudes gives 

(gjjijl] (q 2 \h\2] ^[34] 

(q 2 2) (3A)[(l 2 - k 3 y - ^} ' 1 • ) 

with l 2 = h — k 2 . As for the quadruple cut, it is easily seen that, on this triple cut, 



fell) (, 2 2) ^ (12) 

where we used = I 2 = l\ = fi 2 . The triple-cut integrand then becomes 

[12] [34] // 



(12) (34) [(l 2 -k,Y-ii 



21 



(3.16) 



which, after replacing the three 5^ functions by propagators, integrates to (|3.1j) . where 
we have included an additional (—1) factor following the comments in footnote 6. The 
factor of 2 in (|3.1|) comes from summing over the two "scalar helicities" . The same result 
comes from evaluating the remaining triple cuts. 



We remark that in the case of the quadruple cut we did not even need to insert the 
solutions of the on-shell conditions for the loop momenta into the expression coming from 
the cut. This is not true in general; for example, for the five gluon amplitude discussed 
below the sum over solutions will be essential to obtaining the correct amplitude. 



3.2 The one-loop — h++ amplitude 



The one-loop four gluon scattering amplitude — h++, with a complex scalar running in 
the loop, is given to all orders in e by [38] 



scalar / 1 — 



r,2 + ,3 + ,4" 



2i 



[24? 



st 



(47r) 2 " e [12] (23) (34) [41] u [ su 
t — u . . s — u , , . st T 

-Ms) — —Mt) + 7T- J 4 + k a 



tu 



t 2 



2u 



(3.17) 



We will now show how to derive this result using generalised unitarity cuts. 

First consider the quadruple cut (see Figure 3). The product of tree amplitudes gives 

{l\h\q x \ {q 2 \l 2 \2\ ( 93 |J 3 |3] {q,\k\A} 



i2) (? 3 3) M) 



(3.18) 
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4 + 




3 + 



Figure 3: The quadruple cut for the amplitude 1 2 + 3 + 4 



It is straightforward to show that, on the quadruple cut, 



<? 3 3) (g 4 4) 
(1\IM (q 2 \l 2 \2] 



[23] 



[34] 
(34) 



(l gi ) (g 2 2) [31] 
and hence the quadruple cut in Figure 3 gives 

.2 [3 4] [2 3] 



g(l + ,2+,3 + ,4- 



(3 4) [31] 



(31) 
(23) 

2 (31) 
(2 3) 



[2\h\l) 



[2\h\l) 



In order to compare with (|3.17|) it is useful to notice that 



[34] [23] (31) 



[24] s 



(34) [31] (23) [12] (23) (34) [41] u 
We conclude that the first term in ()3.19|) generates 



Af 



(-, 



[24] s 



(47r) 2 - £ V[ 12 ](23)(34)[41] u 



K 



i ■ 



(3.19) 



(3.20) 



(3.21) 



where the prefactor in ()3.21|) comes from the definition (j3.2j) and ()3.3|) for the function 
K 4 . 

The second term in f!3.19|) corresponds to a linear box integral, which we examine now. 
We notice that the quadruple cut freezes the loop integration on the solution for the cut. 
In the linear box term in ()3.19|) we will then replace l\ in [2|Zi|l) by the solutions of the 
cut, and sum over the different solutions. 
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Specifically, in order to solve for the cut-loop momentum l\ one has to require 

/2 _ ;2 _ ,2 _ ,2 _ 2 
l l — l 2 ~~ 6 3 ~~ 6 4 ~~ A* 5 

Zi = — ki , I2 = Zi — ^2 j Z3 = Z2 — ZC3 , Z4 = Z3 — ZC4 . (3.22) 

In order to solve these conditions, it proves useful [43] to use the four linearly independent 
vectors k\, k 2 , &3 and -ZT, where 



Setting 
one finds 



-K/i : — tfivpa . (3.23) 

Z x = aki + 6/c 2 + ck 3 + dK , (3.24) 



a = h h=X r c = ~h (3 - 25) 



2 ' 



where 



s = {k x + /c 2 ) 2 , t = (k 2 + A; 3 ) 2 , u = (h + fc 3 ) 2 , (3.26) 

and s + t + u = 0. 
Then one has 

[2|li|l> — [21^-^11) = c [2|3|1> = -^[23](31) , (3.27) 

where Z^ denotes the two solutions for the quadruple cut. The square root drops out of the 

calculation (as it should, given that the amplitude is a rational function). We conclude 
that the second term in (j3.19|) gives 7 



f , '5 °A U , (3.28) 



(47r) 2 ~ e V[ 12 ](23)(34)[41] u) 2u 
where 

J n ■= 4[/i 2 ] • (3.29) 
Again, the prefactor in ()3.28|) arises from the definition ()3.3|) . 

In total the quadruple cut (j3.19|) gives 

21 H (K 4 + |- J 4 ) , (3.30) 



(4?r 



i2-e 



where we have again included a factor of two for the contribution of a complex scalar. 
This result matches exactly all the box functions appearing in (|3.17|) . 

7 Recall that in our conventions t := (23) [32]. 
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h 

(a) 



4 + 




(b) 



Figure 4: The two inequivalent triple cuts for the amplitude 1 2 + 3 + 4 + . 



We now move on to consider triple cuts. We start by considering the triple cut in 
Figure 4a, which we label as 2 + , (3 + , 4 + )). It may be shown that this triple cut yields 
the following expression: 



TC(r,2+,(3 + ,4 + )) 



l>' 



[3 4] [2 3] 
(3 4) [31] 



(31) 
(2 3) 



- [2|i!|l> 



1 



-n 



[34] \2\h\l) 
(3 4) (2 3) [31] 



(3.31) 



The first line in (I3.31|) clearly contains the (negative of the) term already studied with 
quadruple cuts - see ()3.19|) (for an explanation of the relative minus sign see footnote 6). 
We now reconsider the linear box term (second term in the first line of (|3.31Jl ). and study 
its Passarino-Veltman (PV) reduction. As we shall see, this box appears also in other 
triple cuts (see ()3.43|) ). 



A 



Let us consider the linear box integral 
d 4 /i d~ 2e fi 



/i 2 It 



(27t) 4 (2tt)-* {il - ^)[(h - k 2 f - ^}[(h -k 2 - k 3 y - n*][(h + k x y - ^} ■ 

(3.32) 

On general grounds the integral is a linear combination of three of the external momenta, 



A» = akt + Pk% + 7 A£ 



(3.33) 
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For the coefficients we find 



a 
P 
7 



(47r) 2 - e 2u 

(4vr) 2 - 2 4 ' 

i 1 
~(47r) 2 - e 2^ 



tJ 4 - 2Ms] 



2J 3 (t) 



(3.34) 



sJ 4 - 2J 3 (s) + 2J 3 (t) 



Taken literally, this means that from the linear box in ()3.31)1 we not only get the J 4 
function but, altogether: 



iM ( st _ i _ , . f _ 
(4tt) z e \2u u u 



(3.35) 



Summarising, the PV reduction of the first line of the triple cut (|3.31|) . lifted to a Feynman 
integral, gives: 

iM 



(4VT) 



2-t 



K 4 + — J 4 - -Ms) + -J 3 (t) 
2u u u 



(3.36) 



The last term in (J3.36|) is clearly spurious - it does not have the right triple cut, and has 
appeared because we lifted the cut-integral to a Feynman integral; hence we will drop it. 
In conclusion, the triple cut (1~,2 + , (3 + ,4 + )) in Figure 4a leads to 



iM 



(47T) 



2-t 



St t 

K 4 + — -h Ms, 

lu u 



(3.37) 



We now consider the last term in (|3.31jh which generates a linear triangle, whose PV 
reduction we consider now. The linear triangle is proportional to 



B^ : = 
On general grounds, 

and hence 



d l\ d jj 



(2tt)4 (2tt)- 2£ ill - ^)[(h - k 2 f - /i 2 ][(/! + hf - /i 2 ] 



B» = 9k1 + rk 



2 ) 



(3.38) 
(3.39) 

\2\B\1) = . (3.40) 

We conclude that the second line in (|3.31j) gives a vanishing contribution, so that the 
content of this triple cut is encoded in (|3.37j) . 

Next we consider the triple cut labelled by ((1 _ , 2 + ), 3 + , 4 + ) and represented in Figure 
46, which gives 



TC((r, 2+),3+,4^ 



[3 4] [2 3] 
(3 4) [3 1] 



(2 3) (2 3) 



1 



[34] (l|31Z 2 -23/ 2 |2] 
(34) (12)[12] (23) [31] 
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(/ 2 + fc 2 ) 2 -/i 2 

(3.41) 



The first term of (j3.41j) clearly corresponds to the function already fixed using quadru- 
ple cuts. The second term can be rewritten as follows. Introducing li := l 2 + k 2 , we have 



(J2) (3|/ 2 |2] = -[2|/ 2 |1> + M [(/ 2 + k 2 f - ^} , (3.42) 



(23) x 1 1 J (23) 
therefore we can rewrite ()3.41|) as 



,[34] [23] / ,,(31) ro| , \ 1 



TC((l-,2 + ),3 + ,4 + ) = - [2\h\l) 



(34) [31] V r (23) L 1 11 7 & + *%,)*- a* 
2 [3 4] / (l|31/ 2 -23/ 2 |2] [23] (31) 



,2 



+ ^ (3 4) V (12) [12] (2 3) [31] [31] (23) J ' (3 ' 43) 

We know already that the PV reduction of the first line of ()3.43|) corresponds to ()3.36|) 
- with the term containing Js(t) removed - so we now study the second line, which will 
give new contributions. 

The second term in the second line corresponds to a scalar triangle, more precisely it 
gives a contribution 

-(SpJiM- < 3 - 44 > 

The first term corresponds to a linear triangle, and now we perform its PV reduction. 
The relevant integral is 

• J (2vr) 4 (2vr)- 2 ^ (Z 2 - /i 2 )[(/ 2 - hf - ^}[(l 2 + h + k 2 f - /i 2 ] " 1 ' ] 
On general grounds, 

C" = A A£ + k (A* + k 2 y . (3.46) 
A quick calculation shows that 



A 



Ms) - -Ms) 



« = 7T^Z7 'Ms) ■ (3.47) 



(47r) 2 ~ e L s W J ' (47r) 2 - e s 

The first term in the second line of ()3.43j) gives then 

b Ja(s) + — j2(s) ) • (3 - 48) 

where A/" is defined in (I3.20J) . Altogether, the second line of ()3.43|) gives 

W ' '1+HV.W + — Ms)) , (3.49) 



(4vr) 2 - e V ^ sJ w s 
whereas from the first line of the same equation we get 



/W f K, + ^J A - -Hs) \ . (3.r,0) 



(4tt) 2 - £ V 2m 4 it 
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where we have dropped the term J 3 (t) for reasons explained earlier. 



We conclude that the function which incorporates all the right cuts in the channels 
considered so far is equal to the sum of ()3.49j) and ()3.5()|) . which gives 

K 4 + — J 4 - - J 3 (s) - 1 + - J 3 (s) + — =- J 2 (s) . (3.51) 



(47T) 2 ~ e V 2u M V s/ ow s 2 

Using — t/w — 1 — u/s = s/u — u/s, f!3.51|) becomes 

iA/* / T r st _ / s it \ T . . u — t r , , 

K 4 + ^-^4 + ) J 3 («) + — 5- J 2(s) . (3.52) 



(47r) 2 ~ e \ 2-u \u s/ s 

To finish the calculation one has to consider the two remaining triple cuts, that is 
(4 + , 1~, (2 + , 3 + )) and ((4 + , 1~), 2 + , 3 + ). These cuts can be obtained from the previously 
considered cuts by exchanging s with t. 

Our conclusion is therefore that the function (including the usual factor of 2) with the 
correct quadruple and triple cuts is: 

2iM ( TJr st _ / s u\ T , . u — t r , s , oro , 
K 4 + — J 4 + J 3 (s) + —5- J 2 (s) (3.53) 



(47r) 2 e \ 2u \u s/ s 

This agrees precisely with ()3.1|) using the identities 

t(w — s) s it s(u — t) t u 

su us 1 tu u t 



(3.54) 



3.3 The one-loop h+ amplitude 

We now turn our attention to the one-loop four point amplitudes with two negative helicity 
gluons. We start by considering the one-loop amplitude *4. 4 calar (l~, 2~, 3 + , 4 + ), which is 
given by [38] 8 

Atree / 7- 1 1 \ 

^scai ar(1 - ; 2 - ; 3+; 4+) = 2 _0_ l^-*-K 4 + -J 2 (t) + -7|- 2e (t)J . (3.55) 

To begin with, we consider the quadruple cut of the amplitude, represented in Figure 5. 
It is given by 

(l|Zi|gi] (2\h\q 2 ] (q 3 \l 3 \3] (g 4 |/ 4 |4] 



[l?i] [2? 2 ] (93 3) (g 4 4) 



(3.56) 



8 Here for simplicity we drop the functions 1\ and -^(O), which are zero in the massless case [38]. We 
also include a factor of two as we are considering complex scalars. 
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4 + 



3 + 



Figure 5: The quadruple cut for the amplitude 1 2 3 + 4 



By choosing q\ = 2, q 2 = 1, qs = 4, g 4 = 3, ()3.56|) can be rewritten as 

t 



s 



where 



A 



tree 



(12) s 



(3.57) 



(3.58) 



(2 3) (3 4) (41) 

Reinstating the four cut propagators and integrating over the loop momentum, (J3.57)) 
gives 

^trce f j. \ 

(3.59) 



^tree / j- 

A ' ~K A 



(4vr 



,2-e 



4 > 



where is defined in fl3.2j) . 



Next we consider triple cuts. We begin our analysis with the triple cut in Figure 6a. 
This yields 



/x 2 [34] (\\h\qy\ [2\h\q 2 ] 4 (12) [3 4] 1 



(3 4)2(Z 2 -3) [l gi ] [2q 2 ] 



[12] (3 4) 2(/ 2 -3) 



(3.60) 



which, upon reinstating the cut propagators and performing the loop momentum integra- 
tion gives 



A*? 00 ft 



(4vr; 



2-e 



(3.61) 



This function had already been detected with the quadruple cut, as discussed earlier. 
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Figure 6: The two inequivalent triple cuts for the amplitude 1 2 3 + 4 + . 

Next we move on to consider the triple cut in Figure 66. This yields 

(1N4] 2 (2\h\ qi ] (ggN3] 
2t(/ 3 -4) [2 qi ] (q 2 3) 

We can re-cast ()3.62j) as follows. Firstly, we write 



t(3 4)[12]2(Z 3 -4) 

(1|4|3](2|1|4]/. 2 
i(34)[12] 

(l|4|3]<2|Z 3 |4]/x 2 
t(34)[12] 



(3.62) 



(jMW! - 2 (1|4|3] 2(/ 3 -4)(l|Z 3 |3] 

Sfe3> - ^ (3 4) (3 4) ' {S - b6) 

and secondly 

(l|W(2jjiM 2 (2|1|4] 2(i 3 -4)(2|l|4] 2(< 3 ■ 4)(2|i 3 |4] 

= " T^r — — m — • (3 - 64) 

The expression 1)3.62)1 becomes a sum of six terms Tj, z = 1, . . . , 6, where 

(1|4|3]<2|1|% 4 
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<2|1|4]<1|* 3 |% : 

La — — 



t(34)[12] 
(2|l|4](l|/ 3 |3]2(/ 3 -4) 



t(3 4>[12] 



= _(l|f,|3](2|f,|4]2(f,.4) 
6 *(3 4)[12] V ; 

Next we replace the delta functions with propagators, and integrate over the loop momen- 
tum. To evaluate the integrals, we use the linear, quadratic and cubic triangle integrals 
in 4 — 2e dimensions listed in the Appendix. The integration of the expressions gives 



/(tree / f 



(4vr) 2 - e V s 



A tree / f 
Atree / f i 

T 3 - -7^z-A-Mt)--Ut) 

(47r) z £ \ s s 

Atree / i 

^ - -(SK -;■*«> 
r « - "(^(- s'*> " (I + 7> r2<(t) " 7^ fc w) ■ < 3 ' 66 > 

We now use (A. 26) in [38] relating J 2 {t) to I 2 (t) and J^ 2e (t), and get 

T 5 + T 6 - -^ ree 0J 2 (O - j/ 2 6 - 2e (*)) . (3.67) 
Adding up the six Tj terms, and including the usual factor of two, we obtain 

r) Atree / f i i \ 

4 ' K,--J 2 {t)--ll-^{t) , (3.68) 



(4tt) 2 - £ Vs s w t 
which precisely agrees with (|3.55j) . 



3.4 The one-loop — I h amplitude 



Now we consider the one-loop amplitude with a complex scalar in the loop, .4.4 Calar (l , 2 + , 3 , 4 + ), 
which is given by [38] 
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^"^(l- 2+3- 4" 



-2 



st 



u z u 



(3.69) 



-2e/ 



st(s - t) 



2^3 



ir 



st(t - s) 



U" 



Ms) + —J 2 {t) + -j J 2 (s) 



tu su 



ts u 



st z 



U" 



u- 




4 + 



Figure 7: The quadruple cut for the amplitude 1 2+3 4 + . 
The relevant quadruple cut is represented in Figure 7, and gives: 

(jjjiN Wgg (3|Z 3 |g 3 ] (^N^] 

(92 2) [3g 3 ] (94 4) 
' , (13)^+(12)(3|/ 1 |2])([24K-[34](3|/ 1 |2] 



[13] (2 4) 



gt/i 4 2s 2 t(|Z 1 |% 2 s 3 t(3|Zi|2] 2 
it 2 + w 2 (3|l|2] + w 2 (3|l|2] 2 



where 



A 



tree 



(13)' 



(3.70) 



(3.71) 



4 ' (12) (23) (34) (41) ' 
Averaging over the two solutions of the quadruple cut we obtain the following expression: 



lA * [ v? * + M 3 ^ + 2^ 



(3.72) 
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After reinstating the four cut propagators and integrating over the loop momentum, ()3.72j) 
gives 



1 



(47T 



JV 



tree 



st 2s 2 t 2 s 3 t 3 
—K^ — J 4 — ——7/4 



u- 



2u 4 



(3.73) 



We now use the identity (A. 26) in [38] ignoring functions that do not have a quadruple 
cut to write this as 

1 AT C ( - /v, + ~—ir At ] . (3.74) 



(47T 




1 + 



Figure 8: The only independent triple cut for the amplitude 1 2 + 3 4 + (the others are 
obtained from this one by cyclic relabeling of the external gluons). 

We now consider triple cuts. There is only one independent triple cut, and we consider, 
for instance, the triple cut in Figure 8, which gives 



(1|/ 3 |4] 2 <3|* 3 |? 2 ] («i|/i|2] 



2t(Z 3 -4) [3q 2 ] ( 9l 2) 



(3.75) 



Using straightforward spinor manipulations, and taking into account properties of the cut 
momenta, one finds that the above expression may be expanded as a product of two sets 
of terms. The first is 



[3 92] 

whereas the second is 

(l|Z 3 |4](gi|Zi|2] 
(<7i 2) 



/i 2 (3|l|4] t(3|/ 3 |4] + 2(/ 3 -4)(3|/ 3 |4] 



[13] 



[13] 



[13] 



/i 2 (l|4|2] , (4|1|2](1|/ 3 |4] 2(/ 3 .4)(l|/ 3 |2] 



(2 4) 



(2 4) 



(2 4) 



(3.76) 



(3.77) 
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The expression ()3.75|) becomes then a sum of nine terms Ri, i = 1, . . . , 9, where 



Ri 



R, 



R* 



Ra 



Re 



Rp, 



R7 



Ra 



(1|4|2](3|1|4]// 
*[13]<2 4)2(Z 3 -4) ' 

(4|l|2](3|l|4](l|Z 3 |4b/ 
t[13](2 4)2(/ 3 -4) 

(3|1|4](1|Z 3 |% 2 
£[13] (2 4) 

(1|4|2](3|/ 3 |4]^ 
[13]<24)2(Z 3 -4) ' 

(4|1|2](3|Z 3 |4](1|Z 3 |4] 
[13](2 4)2(/ 3 -4) 

[13] (2 4) ' 

(l|4|2](3|Z 3 |4b/ 
t[13](2 4) 

(4|1|2](3|Z 3 |4](1|/ 3 |4] 
*[13](24> 



= _(3|Z,|4]<l|Z,|2]2(/,.4) 

t[13](2 4) V ; 

The term R 5 becomes a quadratic box integral when the three delta functions are replaced 
with propagators. We can use the properties of the cut momenta to re-write R$ as a sum 
of terms which will give a box integral, a linear box integral and a linear triangle integral 
as follows, 

(4|l|2][4|31|4]/i 2 t(4|l|2][4|3Z 3 |4] _ (4|1|2][4|3 Z 3 |4] 
5 [13] 2 (2 4)2(Z 3 -4) [13] 2 (2 4)2(Z 3 -4) [13] 2 (24) ' 1 ' 

We now replace the delta functions with propagators and integrate over the cut momenta. 
Note that one must drop any terms without cuts in the t-channel. This must be used for 
all the linear box integrals that appear above. Using the results for the linear box and the 
linear, quadratic and cubic triangle integrals in 4 — 2e dimensions listed in the Appendix 
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gives 



^trcc / g 2 t 2 g 2 t 

^ - -7^f-4^)+4^) 

[4:71 J z € \ vr u z 
^trcc / s 2 t 2 gt 2 

AT e ( st 



fl » - -(^(- 5?'*> + (s? - + s^w) • < 3 ' 80 > 

Now using (A. 26) in [38], and ignoring all terms without cuts in the t-channel, it is easy 
to show that the sum of these nine terms leads to the result 

.4-»'(l-,2+,3-,4+) = -^AT^Kt-^ir' + ^'m (3.81) 

- ^^-Mt) + - 2 m + f it*® + ^/ 2 («)) . 

u 6 u tu u 6 J 

Next, one must also include the corresponding terms coming from the s-channel version 
of the of triple cut in Figure 8. This just yields (j3.81j) with t replaced by s. Combining 
these two expressions, without double- counting the box contributions (which appear in 
both cuts), and including the usual factor of two, one precisely reproduces the amplitude 
for this process ()3.69|) 
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4 The + 



amplitude 



The five-gluon all-plus one loop amplitude, with a scalar in the loop, is given by [50] 



.4 5 (l + ,2 + ,3 + ,4 + ,5+) 



96tt 2 C 5 

where C 5 := (12) (23) (34) (45) (51) and e(abcd) := e^ upa a^b u c p d a 

An expression for the five-gluon amplitude valid to all orders in e appears in [39] 
^ calar (l + ,2 + ,3 + ,4 + ,5 



S12S23 + S23S34 + S34S45 + S45S51 + S51S12 + 4ie(1234) 

(4.1) 

pupa ' 



i e(l-e) 



C 5 (Air) 



2-e 



S23S344 1} ' 8 ^ + S34S454 2) ' 8 ^ 



+ S 4 5S 5 l4 3) ' 8 ^ + S 5 iSi 2 4 4) ' 8 ^ + Si 2 S234 5) ' 8 ^ 



+ 4i(4-2e)e(1234)/; 



10-2e 



The result (|4.1j) is obtained from (|4.2jl by taking the e — > limit, where [39] 



24' 6(1 " eU 



10-2e 
6 



(4.2) 



(4.3) 



1 + 




4 + 



Figure 9: One of the quadruple cuts for the amplitude l + 2 + 3 + 4 + 5 + . 



Here we will find that we can reproduce the full amplitude using only quadruple cuts 
in 4 — 2e dimensions. 
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Let us start by considering the diagram in Figure 9, which represents the quadruple 
cut where gluons 4 and 5 enter the same tree amplitude. The momentum constraints on 
this quadruple cut are given by 

72 _ 72 _ ,2 _ ,2 _ 2 
l l — t 2 — '3 — '4 — l Jj i 

h = U - h , k = h~ k 2 , h = k~ h , U = h~ h~ h ■ (4.4) 

It will prove convenient to solve for the momentum Z 3 , which we expand in the basis of 
vectors ki,k 2 , k 3 and K, where K is defined in ()3.23j) . One finds that the solution of ()4.4|) 
is given by 9 

l 3 = ah + bk 2 + ck 3 + dK , (4.5) 



with 



a = ±- , b = -~ , c = -1 - ± , (4.6) 
2u 2 2u 



St + 4:/J?U 

d = ± v — s^ - ' 

where the kinematical invariants s, t, u are again defined by (j3.26|) . but now s + t + u = 

(h + k 5 ) 2 . 

Considering the diagram in Figure 9, the product of tree-level amplitudes entering the 
quadruple cut can be written as 

Wi Wgl f 45 J (a 7] 

( qi l) ( ft 2> (g 3 3) (45) [(h - ktf - ■ { ■ ] 

Using (|3.10j) . and choosing g 3 = 2, (|4.7J) can be recast as 

4 [12] [45] 1 (2|/ 3 |3] // 4 Tr_(5123/ 3 4) 



(12) (45) (23) (l 3 - k 4 f - (12) (23) (34) (45) (51) (Z 3 - A; 4 )2 - ^ 

// Tr+(123Z 3 43) + Tr+(123Z 3 42) 



(4- 



(12) (23) (34) (45) (51) (/ 3 - k,f - ^ 

Using momentum conservation, and 

Tr + (abed) = 2 (ab)(cd) - (ac)(bd) + (ad) (6c) + ie(abcd) , (4.9) 

it is easy to see that 
Tr + (123/ 3 43) + Tr + (123/ 3 42) = 4(12)(23) _ ^ (34)6(121,3) - (12) e(234Z 3 ) ^ 



9 We notice that, had we solved for li, the solution would have taken the form l|3.24J) with the same 
coefficients a, b, c, d of Ij3.25|l - but with u defined by u = — s — t — (£4 + k^) 2 . 
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We set 



ie(12Z 3 3)(3-4) - ze(234Z 3 )(l • 2) . 



(4.11) 



Now we wish to sum the expression ()4.8|) over the solutions ()4.6|) . including a factor of 1/2. 
Writing these solutions as if = x±y, where y contains the term involving the momentum 
K, it is straightforward to show that 



if 



Tr+(123Z 3 43) + Tr+(123/ 3 42) 
(h ■ h) 



4(1 ■ 2)(2 ■ 3) - 4 



V{x){x-A) - V(y){yA) 



(x ■ 4) 2 - (y ■ 4) 2 



and 



V(x)(x-4) - V{y){y±) 
{x ■ 4) 2 - (y • 4) 2 

Summarising, we have found that 

1 >^ Tr + (123Z 3 43) + Tr+(123Z 3 42) 
2^ (l^h) 

= s 12 s 23 - 4z> 2 e(1234) 



— /i 2 e(1234) 



1 1 

+ 



(If ■ 4) (Is ■ 4) 



(4.12) 
(4.13) 



4(l-2)(2-3) + 2^ 2 e(1234) 
1 



1 



+ 



(h ~ k ±) 2 ~ V 



(tf-4) (/ 3 -4)J 
(4.14) 



From ()4.8j) . we see that the full amplitude in the quadruple cut is obtained by multiplying 
(|4.14J) by —fi 4 /Cs. Next, we lift the cut integral to a full Feynman integral, and get 



-2 



S12 s 23 ~ 4z/i 2 e(1234) 



(/ 3 + - £; 4 ) 2 - M 2 (l^-h) 2 -fi 2 



C 5 (47r) 2 - 
i e(l-e) 

2-e 



r (5),4-2e r 4 



[/i 4 ]s 12 s 23 + 8^/ 5 4 - 2e [/]e(1234) 



C 5 (4tt; 



S12S23-Z, 



(5),8-2e 



+ 4i (4 - 2e) e(1234)J, 



10-2e 



(4.15) 



where the factor of 2 in the first line of ()4.15|) comes from adding, as usual, the two 
possible quadruple cuts of the amplitude (which are equal, since they are obtained one 
from the other by simply flipping all the internal "scalar helicities"). 



Let us now discuss the result we have found. The first term in the last line of ()4.15|) 
gives the Si2S 2 3 term in (|4.2j) . The other quadruple cut diagrams, which come from cyclic 
relabelling of the external legs, will similarly generate the other e(1234)-independent terms 
in ()4.2jl . Finally, the e(1234) term in (|4.15j) - a pentagon integral term - matches the 
e(1234) term in (jOJl . 

Thus we have shown that the five gluon amplitude +++++ may be reconstructed 
directly using quadruple cuts in 4 — 2e dimensions. 
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Appendix A: Tensor Integrals 



In this section we summarise the tensor bubble, tensor triangle and tensor box integrals 
used in this paper. 



The scalar n-point integral functions in D = 4 + 2m — 2e dimensions are defined as 

= = i(-l)<+WP I 1 — i (A.l) 

J (A) L 2 (L-p l ) 2 ---(L-}2i = iPi) 2 

i(-l) n+1 f d 4+2m ld~ 2 ^ 



n 2+m ~* J (P - Pl y -f) •••((/- vi? -/i 2 ) ' 

The higher dimensional integral functions are related to 4 — 2e dimensional integrals with 
a factor fi 2m inserted in the integrand. For m — 1,2 one finds 

J„[/i 2 ] = J n = (-e)7«- 2e , and J n [/i 4 ] ee K n = (-e)(l - e)/*"* . (A.2) 



In our paper we encounter bubble functions with m — 0, 1, triangles with one massive 
external line and m = 0,1, and boxes with four massless external lines and m — 0, 1, 2: 

/,(**) = ^i-pr 1 - , ir-w - - 2e(1 -I) r (1 - 2£) (- f2 )-' • 

(-i)If 2 - = + 0(e) , (-£)(! - ejlt 2 ' = -\ + O(e) . (A.3) 
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Figure 10: Kinematics of the bubble and triangle integral functions studied in this Ap- 
pendix. 

Note that the expressions for the bubbles and triangles are valid to all orders in e, whereas 
for the box functions we have only kept the leading terms which contribute up to 0(e°) 
in the amplitudes. 

We now move on to present the result of the PV reduction for various tensor integrals 
which are relevant for this paper. Note that the expressions are presented in terms of 
scalar n-point integral functions /jf in various dimensions V, specifically in terms of I n , 
In~ 2e and In~ 2e in 4 — 2e, 6 — 2e and 8 — 2e dimensions, respectively. The expressions 
are valid to all orders in e, if I n , In~ 2t and If~ 2e are evaluated to all orders, and the PV 
reductions have been performed in a fashion that naturally leads to coefficients without 
explicit e dependence (the reader may consult [51] for more details on this particular 
variant of PV reductions). 

For the linear and two-tensor bubbles we have (see Figure 10a): 



I 2 [l/g Lg] 





\h + ^/ 2 6 - 2 M (P2 + P3 Y(P2 + P 3 y ■ (A.5) 



(A.4) 
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For the linear, two- and three-tensor triangles (see Figure 106): 

h [14] = -\hp$ + ( - h + , (A.6) 



+ ( - y/I' 26 + I 1 !' 2 ") (^>3 + + ^) • (A.8) 
Finally, for the linear box: 

hM = (±i A -l(i s ( t )-i 3 (s))^-^ 

+ {^h - l(h(t) - I 3 (s))^ , (A.9) 

where, as usual, 1% denote D-dimensional scalar n-point integral functions, s := (P1+P2) 2 , 
t '■= ij>2 + P3) 2 , u := (pi +P3) 2 , and I n is an abbreviation for the (4 — 2e)-dimensional 
integral functions. 
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